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Abstract 

We study the effects of color superconductivity on the structure and formation of 
compact stars. We show that it is possible to satisfy most of recent observational 
boundaries on masses and radii if a diquark condensate forms in a hybrid or a quark 
star. Moreover, we find that a huge amount of energy, of the order of 10^^ erg, can 
be released in the conversion from a (metastable) hadronic star into a (stable) 
hybrid or quark star, if the presence of a color superconducting phase is taken into 
account. Accordingly to the scenario proposed in Astrophys. J. 586(2003) 1250, the 
energy released in this conversion can power a Gamma Ray Burst. This mechanism 
can explain the recent observations indicating a delay, of the order of days or years, 
between a few Supernova explosions and the subsequent Gamma Ray Burst. 
PACS: 26.60.+C, 26.50.+X, 12.38.Mh, 97.60.Jd 
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1 Introduction 



The new accumulating data from X-ray satellites provide important informa- 
tion on the structure and formation of compact stellar objects. Concerning the 
structure, the new data fix rather stringent constraints on the mass and the 
radius of a compact star. These data are at first sight difficult to interpret in a 
unique and self-consistent theoretical scenario, since some of the observations 
are indicating rather small radii and other observations are indicating large 
values for the mass of the star. 

Concerning the formation scenario, crucial information are provided by the 
very recent observations of Gamma- Ray Bursts (GRB), indicating the possi- 
bility that some of the GRBs are associated with a previous Supernova ex- 
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plosion, with a delay between the first and the second explosion of the order 
of days or years [1,2]. These observations could be explained associating the 
second explosion with the conversion of a (metastable) hadronic star (HS) into 
a more stable stellar object made at least in part of deconfined quark matter 
(QM). In the scenario proposed in Ref.[3] (see also references therein), the HS 
can be metastable due to the presence of a non-vanishing surface tension at 
the interface separating hadronic matter (HM) from QM. The nuclcation time 
(i.e. the time to form a critical-size drop of quark matter) can be extremely 
long if the mass of the star is small. Via mass accretion the nucleation time 
can be dramatically reduced and the star is finally converted into the stable 
configuration. 

In recent years, many theoretical works have investigated the possible forma- 
tion of a diquark condensate in quark matter, at densities reachable in the 
core of a compact star [4,5,6]. The formation of this condensate can deeply 
modify the structure of the star [7,8,9]. 

In this Letter we show that it is possible to satisfy the existing boundaries on 
mass and radius of a compact stellar object if a diquark condensate forms in a 
Hybrid Star (HyS) or a Quark Star (QS). Moreover, the formation of diquark 
condensate can significantly increase the energy released in the conversion 
from a purely HS into a more stable star containing deconfined QM. 



2 Equation of state of beta-stable matter 

The EOS appropriate to the description of a compact star has to satisfy beta- 
stability conditions. Moreover two charges are conserved, the baryonic and 
the electric one. These conditions need to be satisfied in the hadronic, in the 
quark and in the mixed phase. Although electric charge neutrality and beta- 
stability are easy to impose for non-interacting quark matter, these conditions 
are highly non-trivial when a diquark condensate can develop. 

Concerning the hadronic phase we use the relativistic non-linear Glendenning- 

Moszkowski model (GM1-GM3) [10]. At very low density we have used the 
Ncgclc-Vauthcrin [11] and the Baym-Pethick-Sutherland [12] EOS. For the 
quark matter phase we adopt a MIT-bag like model in which the formation 
of a diquark condensate is taken into account in a simple and effective way 
which will be described below. To connect the two phases of our EOS, we 
impose Gibbs equilibrium conditions. When the Gibbs conditions are applied 
in presence of more than one conserved charge, the technique developed by 
Glcndenning has to be adopted [13,14] and the pressure need not to be con- 
stant in the mixed phase. Therefore a finite volume of the star can be occupied 
by the mixed phase, what is crucial for the stability of the star. 
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It is widely accepted that the Color-Flavor Locking phase (CFL) is the real 
ground state of QCD at asymptoticly large densities but, at the scales which 
are involved in a compact star, there is yet uncertainty about the presence 
of this phase and in particular about the transition from superconducting 
QM to HM. There are different possible scenarios depending on the value 
of the strange quark mass. A direct transition from CFL to HM is possible 
for a small value of m^. Alternatively, an intermediate window of Crystalline 
Color Superconductivity phase can develop [15,16] or a phase can be formed 
in which the Cooper pairing involves only quarks of the same flavor [17]. 
A pure two flavor color superconductivity (2SC) phase is ruled out due to 
the high free energy cost resulting from the requirement of color and electric 
neutrality [17,18]. A mixed 2SC-CFL phase could exist but it seems to be 
unstable due to Coulomb and surface effects [19]. In our model the CFL phase 
is connected directly to HM through a flrst order phase transition, what seems 
to be consistent with the use of a small value for m^, of the order of 150 MeV. 

In this Letter we are interested in the bulk properties of a compact star. We use 
therefore a schematic model which takes into account in a simple and effective 
way the main characteristics of the EOS of quark matter in the presence of 
diquark condensation. The main aim of our work is to discuss the dependence 
of the structure of the star on the numerical value of three crucial parameters, 
namely the height and position of the maximum of the diqiiark gap and the 
value of the pressure of the vacuum B of the MIT bag model. 

We adopt the scheme proposed in Refs.[7,20] where the thermodynamic po- 
tential is given by the sum of two contributions. The first term corresponds 
to a "fictional" state of unpaired quark matter in which all quarks have a 
common Fermi momentum chosen to minimize the thermodynamic potential. 
The other term is the binding energy of the diquark condensate expanded up 
to order (A//x)^. In Ref.[7] the gap is assumed to be constant, independent 
on the chemical potential. In the present calculation we consider a ji depen- 
dent gap resulting from the solution of the gap equation. We describe the 
superconducting phase using the first work of Afford, Rajagopal and Wilczek 
in which the Color-Flavor Locking phase was introduced [4] by considering 
three massless flavors. This approximation is a sensible one as long as is 
small in comparison with the quark chemical potential [15]. The quark-quark 
interaction is described by a NJL-like Lagrangian with an effective coupling 
constant K and a form factor which mimics the asymptotic freedom of QCD. 



parameters of this model. 

The CFL phase is characterized by the existence of two order parameters 
(singlet) and Ao (octet) which are the solutions of two coupled gap equations. 
The value of the coupling K is fixed, as in Ref.[4], by imposing that the NJL 
Lagrangian with the same K gives, through chiral symmetry breaking at zero 




where w and A are free 
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chemical potential, a reasonable value for the chiral gap A^. This procedure 
in turn fixes the maximum of the superconducting gap. In this work we use 
values for K giving a ranging from 300 to 600 MeV and a corresponding 

maximum of the "effective gap" A(/x) = F^(/^)^[8(Ao(/x))2 + (A,(/i))2]/12 
varying from 70 to 150 MeV. This definition of A corresponds to the gap used 
in Ref.[7] if A^ = 2|Ao|, as assumed in that work. In Fig.l we display the 
"effective gap" as a function of the chemical potential. In our calculation we 
have not really solved self-consistently the coupled equations for the chiral 
and the superconducting gap. The result of microscopic calculations, like the 
ones of Refs. [18,19] indicates that the two gaps are mutually exclusive. In 
particular, the superconducting gap is suppressed at low /i. Therefore we do 
not consider realistic the parameters corresponding to gaps Ai and A2, and 
they are mainly introduced to illustrate the effects connected with the position 
IJ"max of the maximum of the gap. 

In our model, confinement is schematically described by introducing the MIT 
bag constant B. Moreover, the pressure and the energy density are modified by 
the contributions of the electrons, which are necessary in the mixed phase ^ : 



P = -QcFlifi) -B- 1]^'^^*^°"^(/Xe) (1) 

E/V = OcFl(/x) + A^P + 5 + n^^^^^'-°"^(/Xe) + ^JiePe , (2) 

where 

^cfM ^-J e{k-ix)dk + ^ e(Jk^ + ml -i,)dk- (3) 



with 



mi 



(4) 



and the quark density p is calculated numerically by deriving the thermody- 
namic potential respect to 

In Fig. 2 we show the EOS with and without color superconductivity. The effect 
of the gap is to increase the pressure of paired QM respect to the unpaired 
QM at a fixed chemical potential. This extra pressure reduces the values of the 
critical densities. Comparing the curves for A = and A 7^ we can see that 
the CFL EOS is softer than the unpaired QM EOS at low density and stiffer 
at high density. This will have important consequences in the M-R curves. 

^ In the pure CFL phase the contribution of electrons vanishes due to the electrical 
neutrality enforced by the existence of the gap [21]. 
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Fig. 2. Pressure versus baryonic density. HM indicates a purely hadronie EOS, MP 
a mixed-phase of hadrons and quarks and QM pure quark matter. The effect of a 
non-vanishing superconducting gap is displayed. 



3 Masses and radii of compact stellar objects 



In Fig. 3 we have collected most of the analysis of data from X-ray satellites, 
concerning masses and radii of compact stellar objects [22,23,24,25,26,27,28]. 
Most of the data have been obtained very recently due to the technological 
progresses in the field of X-ray detectors. Although some (or all) of the data 
analysis are controversial, since they depend on specific assumptions on the 
structure of the X-ray source, we do think that these observational results 
deserve to be carefully discussed. 
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Fig. 3. Mass-radius plane with observational limits and a few representative theo- 
retical curves: thick solid line indicates CFL quark stars, thick dot-dashed line CFL 
hybrid stars, thick-dashed line hadronic stars (see text). The observational limits 
come from: (a) Sanwal ct al. 2002 [22], (b) Cottam et al. 2002 [23], (c) Quaintrell 
et al. 2003 [24], (d) Heinke et al. 2003 [25], (e),(g) Dey et al. 1998 [26], (f) Li et al. 
1999 [27], (h) Burwitz et al. 2002 [28]. 



Observing Fig.3, we notice that the constraints coming from a few data sets 
(labeled "e" , "f " ^ "g" and maybe also constraint "h" ^ ) indicate rather un- 
ambiguously the existence of very compact stellar objects, having a radius 
smaller than ~ 10 km. At the contrary, at least in one case ("a" in the figure), 
the analysis of the data suggests the existence of stellar objects having radii 
of the order of 12 km or larger, if their mass is of the order of 1.4 Mq. In this 
analysis one has also to take into account that it is difficult from an astro- 
physical viewpoint to generate compact stellar objects having a mass of the 
order of one solar mass or smaller. Therefore the most likely interpretation of 
constraint "a" is that the corresponding stellar object does not belong to the 

^ A very recent reanalysis of the data of the pulsar SAX J1808. 4-3658, discussed in 
Ref. [27], seems to indicate slightly larger radii, of the order of 9-10 km for a star 
having a mass of 1.4-1.5 Mq [29]. 

^ The data at the origin of constraint "h" have been discussed in many recent 
papers. In Ref. [30] an indication for an even more compact stellar object can be 
found. Anyway, the so-called thermal radius obtained in these analysis could be 
significantly smaller than the total radius of the star. 
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same class of objects which have a radius smaller than ~ 10 km. Concerning 
constraint "b", its interpretation is less clear, since it can be satisfied both 
with a very compact star or with a star having a larger radius. The apparent 
contradiction between the constraints "e", "f", "g" and the constraint "a" 
can be easily accommodated in our scheme, since it can be the signal of the 
existence of metastable purely hadronic stars which can collapse into a stable 
configuration when deconfined quark matter forms inside the star. In the next 
Section we will discuss the possible relation between this transition and at 
least some GRBs. 

Finally, constraints ("c" and "d") do not provide stringent limits on the radius 
of the star, but they put strong constraints on the lower value of its mass. 
Constraints "c" and "d" are very important, since it is in general not easy 
to obtain solutions of the Tolman-Oppenheimer-Volkoff equation having both 
large masses and very small radii. As we will see, the existence of an energy gap 
associated with the diquark condensate helps in circumventing this difficulty, 
since the effect of the gap is to increase the maximum mass of QSs or of HySs 
having a huge content of pure quark matter, as shown in Fig.4. 

In the upper panel of Fig.4 we show that a CFL HyS has a smaller radius 
respect to an unpaired QM HyS, if the mass of the star is smaller than ~ 1.35 
Mq. This can be explained observing the EOSs shown in Fig. 2. For a low mass 
star, the central density lies in the region of the P — p plane where the CFL 
EOS is softer than the unpaired QM EOS and therefore the radius of the star 
is smaller. For a large mass star, on the contrary, the central density is in the 
region where the CFL EOS is stiff er than unpaired QM EOS and therefore 
the radius of the star is larger. It is also worth remarking that if the value of 
the gap is increased the amount of QM in the star also increases. Therefore, 
for large values of the gap, heavy HySs have a shape more and more similar 
to the shape of pure QSs, which are finally obtained by a further increase of 
the value of the gap^ (see also the lower panel of Fig.4). All these results are 
totally consistent with the ones obtained in Ref.[7]. 

Concerning the value of the chemical potential Umax-, which corresponds to 
the maximum of the gap, we get relevant modifications to the mass-radius 
relation if 0.3 GeV < Pmax ^ 0.6 GeV. For larger values of Pmax the effect of 
the gap is negligible. Low values of Pmax, of the order of 0.3 GeV, correspond 
to interesting stellar configurations, but are difficult to justify at the light of 
results hke the ones presented in Ref.[19] (see discussion at the end of Sec. 2). 
In the upper panel of Fig.4 we display results for three values of Pmax, whose 
corresponding gaps are shown in Fig.l. We also show that, using B^^^ = 170 
MeV, only a big value of the gap, located at a not too large density (A4) 



^ Notice that the existence of a large CFL gap is not constrained by the traditional 
argument concerning the stability of Fe against decay into two-flavor QM. 



7 



allows the formation of a QS, while for the other gaps HySs are obtained. 

In Fig. 3 we show a few theoretical M-R relations which correspond to the 
scenario we are proposing. More precisely, we show a thick-dashed line corre- 
sponding to HSs (GMl), a thick dot-dashed line corresponding to HySs (GMl, 
B^/^ = 170 MeV, A2) and a thick solid line corresponding to QSs (GMl, 
j^i/i _ -^^Q MeV, A4). Both the HyS and the QS lines can satisfy essentially 
all the constrains derived from observations (concerning the constraint "f" see 
footnote (2)). Concerning the constraint "a", it is probably better satisfied by 
the HS line than by the HyS or QS hnes, which would give stars having a mass 
smaller than ~ 1.2Mq. In conclusion, in our scheme most of the compact stars 
are either HySs or QSs having a mass in the range 1.2 — I.SMq and a radius 
~ 8.5 — 10 km. Metastable HS can exist. As we will see in the next section 
their mass is probably smaller than ~ l.SM©. 



4 Nucleation time and energy released 

In the model we are discussing the formation of QSs or HySs is due to the 
conversion of a purely HS into a more compact star in which deconfined QM is 
present. An HS can be metastable if a non- vanishing surface tension is present 
at the interface between HM and QM. The process of quark deconfinement 
can be a powerful source for GRBs and it can also explain the delay between a 
supernova explosion and the subsequent GRB observed in a few cases [1,2] ^ . In 
the scenario proposed in Ref.[3], the central density of a pure HS increases, due 
to spin down or mass accretion, until its value approaches the deconfinement 
critical density. At this point a spherical virtual drop of QM can form. The 
potential energy for fluctuations of the drop radius R has the form [33] : 

4 

U{R) = -TTR^nqiiig - Hh) + 47r(Ti?2 + M-fR (5) 

where Uq is the quark baryon density, jih and are the hadronic and quark 
chemical potentials, all computed at a fixed pressure P, and a is the surface 

^ A possible mechanism explaining these GRBs is the supranova model [31]. In this 
model, the GRB is the result of the collapse to a black hole of a supramassive fast 
rotating NS, as it loses angular momentum. According to this model NS is pro- 
duced in the Supernova explosion preceding the GRB event. The initial baryonic 
mass of the NS is assumed to be above the maximum baryonic mass for non-rotating 
configurations. As noticed in Ref. [32] in these collapse too much baryonic mate- 
rial is ejected and thus the energy output is expected to be too small to produce 
GRBs. Moreover, the supranova model seems to produce GRBs which are too short 
compared with the observed durations. 



8 




M/Mo 

2.5 r 




Fig. 4. A few theoretical mass-radius relations are shown. HS indicates purely 
hadronic stars. The other lines correspond either to hybrid or quark stars, depending 
on the value of the gap for a given value of B. 



tension for the interface separating quarks from hadrons. Finally, the term 
containing 7 is the so called curvature energy. For a we use standard values 

from 10 to 40 McV/fm^ and we assume that it takes into account, in a effective 
way, also the curvature energy. The value of a was estimated in Ref.[34] to 
be ~ 10 MeV/fm^. Values for a larger than ~ 30 MeV/fm^ are probably not 
useful at the light of the result of Refs. [35,36]. 
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To compute the time needed to form a bubble of quarks having a radius larger 
than the critical one, we use the technique of quantum tunnehng nucleation. 
We can assume that the temperature has no effect in our scheme: for values 
of S^/^ ~ 160 — 180 MeV, which we use in this Letter, the critical density 
pi separating pure HM from mixed phase is larger than 4po for ZjA ~ 0.3, 
i.e. for an isospin fraction typical of a newly formed and hot proto-neutron 
star [37]. This critical density typically exceeds the central density of hot and 
not too massive stars. Therefore the mixed phase can form only when the 
star has deleptonized and its temperature has dropped down to a few MeV 
[38] . When the temperature is so low, only quantum tunneling is a practicable 
mechanism. The calculation proceed in the usual way: after the computation 
(in WKB approximation) of the ground state energy Eq and of the oscillation 
frequency of the virtual QM drop in the potential well U{R), it is possible 
to calculate in a relativistic frame the probability of tunneling as 



Po = exp[ —] (6) 



where 

R+ 



A{E) = 2 J dR^[2M{R) + E- U{R)][U{R) - E] . (7) 



R- 



Here R± are the classical turning points and 

2 



M{R) = 47rpJl-'^yR' , (8) 



Ph being the hadronic energy density and n^, riq are the baryonic densities at 
a same and given pressure in the hadronic and quark phase, respectively. The 
nucleation time is then equal to 

T = {voPoN,)-^ , (9) 



where Nc is the number of centers of droplet formation in the star, and it 
is of the order of 10''^ [39]. In the calculation of nucleation times we neglect 
the effects of color superconductivity. We assume that the CFL gap cannot 
form until the radius of the quark drop has increased enough and therefore 
the energy released in the pairing process has not to be taken into account 
when computing the nucleation time. A support to our assumption can be 
found in Ref.[40], where the authors have investigated finite size effects on 
the formation of a 2SC gap, projecting onto states of defined baryon number 
and onto color singlets. If the radius of the quark nugget is smaller than a 
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critical length (of the order of 1.5-2 fm in their case) the magnitude of gap is 
drastically reduced. A similar calculation for the CFL phase has not yet been 
done, but it is reasonable to assume that due to the "locking" between color 
and flavor in the CFL phase the color projection will yield even larger effects. 

Let us recall once again the astrophysical scenario we have in mind. In a few 
cases a delay of the order of days or years between the Supernova explosion and 
the subsequent GRB have been postulated to explain the astrophysical data 
on the GRBs. In the scheme we are discussing, this delay is due to the forma- 
tion of a metastable HS having a relatively small mass. The nucleation time, 
computed using Eq. (9), can be extremely long if the mass of the metastable 
star is small enough. Via mass accretion the nucleation time can be reduced 
from values of the order of the age of the universe down to a value of the 
order of days or years. We can therefore determine the critical mass Mcr of 
the metastable HS for which the nucleation time corresponds to a fixed small 
value (1 year in Table 1). 

In Table 1 we show the value of Mcr for various sets of model parameters. In the 
conversion process from a metastable HS into an HyS or a QS a huge amount 
of energy AE is released. AE is the difference between the gravitational mass 
of the metastable HS and that of the final HyS or QS having the same baryonic 
mass. We see in the Table that the formation of a CFL phase allows to obtain 
values for AE which are one order of magnitude larger than the corresponding 
AE of the unpaired QM case (A = 0). Moreover, we can observe that AE 
depends both on magnitude and position of the gap. 

Finally let us comment on the dependence of the results of Table 1 on the value 
oi B. As we can see, if B^^^ < 160 MeV the value of Mcr is very small and it is 
unhkely that a metastable HS having a mass M < Mcr can be obtained from 
a Supernova explosion. In that case all compact stars would be either HySs or 
QSs. If, at the contrary, B^^^ > 170 MeV, the value of M^. is so large that a 
compact star can be only a (metastable) HS since, after the conversion from 
HM into QM, the HS collapses into a Black Hole (BH). Therefore, only if 160 
MeV < < 170 MeV a GRB can be generated within the mechanism we 
are proposing. 



5 Conclusions 

We have studied the effect of color superconductivity on the EOS of quark 
matter and on the mass-radius relation for hybrid and quark stars. Compar- 
ing the theoretical curves with recent analysis of observational data, we find 
that color superconductivity is a crucial ingredient in order to satisfy all the 
constraints coming from observations. The most difficult problem posed by 
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the astrophysical data is the indication of the existence of stars which are 
both very compact (R < 9-10 km) and rather massive (M > 1.7 Mq). We can 
satisfy these constraints either with hybrid or quark stars. In particular, con- 
cerning hybrid stars, the gap increases significantly the maximum mass of the 
stable configuration, while keeping the corresponding radius < 10 km. These 
findings are in agreement with the ones of Ref.[7], where a chemical potential 
independent gap was used. Anyway, we have shown in this Letter that, to ob- 
tain very large masses for hybrid stars, the maximum of the superconducting 
gap need to be located at a value of the chemical potential which is probably 
too small. 



The superconducting gap affects also deeply the energy released in the con- 
version from hadronic star into hybrid or quark star. To explain recent obser- 
vations indicating a delay between a Supernova explosion and the subsequent 
Gamma Ray Burst [1,2], in Ref.[3] it has been proposed to associate the sec- 
ond explosion with the transition from a metastable hadronic star to a stable 
star containing deconfined quark matter. In this Letter we have shown that 
the energy released, which will power the Gamma Ray Burst, is significantly 
increased by the effect of the superconducting gap and it can reach a value of 
the order of 10^^ erg. 

To satisfy the constraints on masses and radii of the compact stellar objects 
and to obtain a huge energy from the conversion of the metastable hadronic 
star into a quark or hybrid star, rather stringent limits on the parameter 
values have to be imposed. More explicitly, the pressure of the vacuum has 
to be in the range 160 MeV < B^^^ ^170 MeV and the value of the gap 
has to have a maximum A ~ 0.15 GeV at a chemical potential 0.4 GeV 
^ pi-max ^ 0.6 GeV. All these parameters are compatible with the results both 
of hadronic physics calculations and of microscopic studies of superconducting 
quark matter. Finally, let us remark that for values of B in the indicated range, 
hybrid stars are obtained for unpaired quark matter, while, in most cases, the 
formation of the superconducting gap yields quark star configurations. 

Very recently, a new candidate for a delayed Supernova - Gamma Ray Burst 
association, with an estimated delay of the order of months, has been proposed 
in Ref.[41]. If confirmed, this observation would constitute an important sup- 
port to models like the one we have suggested. 



It is a pleasure to thank I. Bombaci, G. Fiorentini, F. Frontera, D. Guetta 
and S. Zane for many useful discussions. 
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Table 1 

Energy released A.E (measured in foe=10^^ erg) in the conversion to hybrid or quark 
star (labeled with a •), for various sets of model parameters, assuming the hadronic 
star mean life-time r = 1 yr (see text) . Mcr is the gravitational mass of the hadronic 
star at which the transition takes place, for fixed values of the surface tension a and 
of the mean life-time r. BH indicates that the hadronic star collapses to a Black 
Hole. We indicate with a dash (-) situations in which the Gibbs construction does 
not provide a mechanically stable EOS. 
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